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ABSTRACT 

The properties of both observed galaxies and dark matter haloes in simulations depend on their 
environment. The term "environment" has, however, been used to describe a wide variety of 
measures that may or may not correlate with each other Popular measures of environment 
include, for example, the distance to the A^* nearest neighbour, the number density of objects 
within some distance, or, for the case of galaxies only, the mass of the host dark matter halo. 
Here we use results from the Millennium simulation and a semi-analytic model for galaxy for- 
mation to quantify the relations between different measures of environment and halo mass. We 
show that the environmental parameters used in the observational literature are in effect mea- 
sures of halo mass, even if they are measured for a fixed stellar mass. The strongest correlation 
between environmental density and halo mass arises when the number of objects is counted 
out to a distance of 1 .5 - 2 times the virial radius of the host halo and when the galaxies/haloes 
are required to be relatively bright/massive. For observational studies this virial radius is not 
easily determined, but the number of neighbours out to 1 - 2 /z^'Mpc gives a similarly strong 
correlation with halo mass. For the distance to the A^* nearest neighbour the (anti-)correlation 
with halo mass is nearly as strong provided N ^ 2. We demonstrate that this environmental 
parameter becomes insensitive to halo mass if it is constructed from dimensionless quantities. 
This can be achieved by scaling the minimum luminosity/mass of neighbours to that of the 
object that the environment is determined for and by dividing the distance to a length scale 
associated with either the neighbour or the galaxy under consideration. We show how such 
a halo mass independent environmental parameter can be defined for both observational and 
numerical studies. The results presented here will help future studies to disentangle the effects 
of halo mass and external environment on the properties of galaxies and dark matter haloes. 
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1 INTRODUCTION 

The formation and evolution of galaxies depends on both internal 
and external processes ('nature vs. nurture'). Among the internal 
processes are radiative cooling, the formation of and feedback from 
stars, and accretion of gas onto and feedback from super-massive 
black holes. It is generally assumed that halo mass is the fundamen- 
tal parameter that drives the internal processes for isolated galax- 
ies. External processes are important because galaxies do not live 
alone in the Universe. Galaxy interactions can induce gravitational 
torques that can significantly alter the angular momentum structure 
of the matter in galaxies. This can for example lead to a starburst 
or to more rapid accretion onto the central black hole, which may 
trigger a quasar phase. Smaller galaxies may accrete onto the halo 
of a more massive galaxy. As a galaxy moves through the gaseous 
halo of a more massive galaxy, it may lose gas due to ram pressure 
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forces. Winds and radiation from nearby neighbours may also affect 
the evolution of a galaxy. It is still an open question to what extent 
the properties of galaxies are determined by internal and external 
processes. 

Even if halo mass were the only driver of galaxy evolution, 
galaxy properties would still be correlated with their environment. 
Because peaks in the initial Gaussian density field cluster together, 
more massive galaxies will live close to each other ('galaxy bias'). 
A correlation between surrounding galaxy density and internal 
galaxy properties therefore does not necessarily imply a causal re- 
lation between the two. 

Early, analytic models pr edicted that the clustering of haloes 

lCole&Kaiseilll989l : 
19991) found : 



depends only on t h eir mass ( Kaisej 1984 

1 k „ , semi- 
analytic models of galaxy formation that, to first order, the only 
property of a dark matter halo that correlates with the (pro- 
jected) number density of surrounding galaxies is halo mass. 
Other properties like spin parameter, formation time and con- 
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centration do not strongly depend on the surrounding dark mat- 
ter density. However, later papers have shown that clustering 
also depends on properties like formation time, concentration, 
substructure content, spin and shape, even for fixed mass (e.g. 

1 1 I ' 1 ■ r-" 1 ■ I 



Gao et al.'2005';'Harker et al.'2006';"Wechsler et al.'2006VBett etdj 
2007; Gao & White 2007; Jing et al. 2007; Maccio et al. 2007F 
Wetzel et alj |2007| ; lAngulo et alJ l2008l ; iFaltenbacher & White! 



2OIOI) . All dependencies other than the one with halo mass are. 



however, second-order effects. The formation time and the halo 
merger rate are also fou nd to depend only weakl y on environment 
jGottlober et al.l 
l2009l ; lHahn etal 



2001 



[Sheth & Tormeni |2004 iFakhouri & Mai 

j |2009h . Th ese findings are con sistent with the hy- 

drodynamical simulations of lCrain et al.l ( l2009h . who found that all 
the variations in the properties of simulated galaxies with environ- 
ment can be accounted for by the dependence of the halo mass 
function on environment. 

For both observations and simulations it is difficult to disen- 
tangle the effects of halo mass from those of the external environ- 
ment. The two are correlated (higher mass haloes live, on average, 
in denser environments) and finding an environmental parameter 
that does not correlate with halo mass is non-trivial. Of course, the 
mass of the dark matter halo hosting a galaxy is important for the 
evolution of that galaxy, so halo mass is as good an environmen- 
tal parameter as any other. One would, however, like to be able to 
distinguish halo mass (the "internal environment") from the envi- 
ronment on large scales (the "external environment"). It is not a 
priori clear whether the environmental parameters used in the liter- 
ature measure halo mass, and if so, whether they measure only halo 
mass, or whether they are also, or predominantly, sensitive to the 
external environment. 

Observationally, halo mass is hard to determine. Group cat- 
alogues, abundance matching, clustering, and gravitational lensing 
all provide statistical measures of halo mass. Nonetheless, most ob- 
servational data sets will have to do without dark matter halo mass 
and define environmental parameters based on the distribution of 
visible matter (usually stellar luminosity). 

Many observational studies have, nevertheless, investigated 
the effect of the environment on the physical properties of galax- 
ies. In general, in higher densit y environm ents galax ies form 



their stars earlier and fast er (e.g. [ Lewis etal.l 2002: Bal drv et al 
2004'; 'Balogh et al."2004a''5; 'Kauff mann et al.ll2004 ; iThomas et al 



2005; Smith etal. 2006) and galaxy morphologies become more 
(pressure support dominated) early t ype, a s opposed to (rotation 
dominated) late typ e (e.g. IDressleif[l980l; lOressler et"^ Il997l ; 
IWilman etal1l2009h . From observations alone it is very hard to 
judge whether these trends are driven mostly by halo mass or 
whether other halo properties and/or large-scale environment play 
an important role. As in observations environment is usually con- 
trasted with stellar mass (rather than halo mass), such observation- 
ally based distinctions between stellar mass and environment may 
tell us more about the stellar mass - halo mass relation than about 
the difference between external environment and halo mass. 

In simulations, halo mass (and other halo parameters) are 
readily available. From simulations much 'cleaner' definitions of 
environment can be obtained, as the distance between objects is 
known in three dimensions, contrary to observations which can 
only provide precise separations perpendicular to the line of sight. 
Radial velocity differences give an indication of separations along 
the line of sight, but peculiar velocities complicate their interpreta- 
tion. 

Many different measures of environment have been used in 
the literature. Some are closely related by construction, while the 



relation between others is more obscure. In this paper we com- 
pare several popular indicators of environments. The aim is to in- 
vestigate which indicators correlate strongly with each other and 
with halo mass and which ones do not. We measure environmen- 
tal parameters using a semi-analytic model for galaxy formation 
(De Lucia & Blaizot 2007) constructed on the merg er tree of dark 
matte r haloes formed in the Millennium Simulation dSpringel et al.l 
l2005l) . so that we also have halo masses available. We will present 
environmental parameters that measure halo mass, but are insensi- 
tive to external environment, along with environmental parameters 
that are insensitive to halo mass. These can be used for studies that 
aim to separate the effect of halo mass and external environment. 
We will show that most of the environmental indicators used in lit- 
erature are in effect measures of halo mass. In the remainder of the 
paper we will use the term 'environment' whenever we quantify 
distances to nearby galaxies, surrounding galaxy densities etc., but 
never when referring to halo mass, in order to clearly distinguish 
the two. 

This paper is organised as follows. Section |2] gives a short 
overview of the literature on environmental parameters, both from 
observations and simulations. In Section[3]we investigate how often 
used environmental parameters correlate with host halo mass. The 
strength of the correlation with halo mass depends on the distance 
scale used in the environmental parameters, as we will show in Sec- 
tion |4] In Section [5] we discuss how to construct an environmental 
parameter that is insensitive to halo mass. Finally, we conclude in 
Section [6] In appendix |A] we provide fitting functions for the host 
halo mass as a function of galaxy environmental parameters. 



2 POPULAR ENVIRONMENTAL PARAMETERS 

The study of the effect of environment on the evolution of galaxies 
has undergone considerable progress through large ga laxy surveys 
like the Sloan Di gital Sky Survey (SPSS; Stoughto n et al.ll200% 
and (z)COSMOS [scoville et alj2007l ; lLillv et alj2007 !). Many dif- 
ferent definitions of environmental density exist. Observationally, 
the density around galaxies must usually be based on the distribu- 
tion of the galaxies themselves, as the distribution of mass is very 
hard to measure. Two slightly different measures are used very of- 
ten: the number of galaxies within a fixed distance and the distance 
to the A'* nearest neighbour. Table [T|contains a short summary of 
the literature on the environmental dependence of galaxy proper- 
ties, both from observations and from simulations. We will expand 
on these in this section in the next section we will study some of 
these in more detail using the galaxy catalogues from the Millen- 
nium database. 

For the environmental parameters it is important, as we will 
show below, whether the masses of the other galaxies used to mea- 
sure the environment have a fixed physical lower limit (or luminos- 
ity), or whether the minimum mass is a fixed fraction of the mass 
of the galaxy one wants to know the environment of. It also matters 
whether the distance out to which the environment is measured is 
fixed in absolute terms or whether it is fixed relative to some length 
scale related to the galaxy in question (e.g. the virial radius of its 
host halo). In Table[T]we indicate for each environmental parameter 
listed (described in the first column) out to what distance (or dis- 
tance equivalent parameter) the environment is measured (second 
column), and whether the minimum mass/luminosity of the galax- 
ies used for the environmental estimate is fixed in absolute terms or 
whether it is a fixed fraction of the mass/luminosity of the galaxy 
in question (if applicable, third column). The last column lists ref- 
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erences to papers employing the parameter. From Table[T]it is clear 
that very few papers take minimum masses of neighbours and/or 
distances relative to properties of the galaxy's host halo. 

Two main classes of observational parameters can be identi- 
fied: those which measure the number of galaxies out to a given 
distance, and those that measure the distance out to a given A'"' 
neighbour. Note that the number of galaxies out to a distance r 
is equivalent to the local number density of that same sample of 
galaxies smoothed on the scale r with a top-hat filter. Similarly, the 
distance to the nearest neighbour, is equivalent to the local 
number density of galaxies smoothed on scale In higher density 
regions the Af"" nearest neighbour is, on average, closer by and the 
scale on which the environment is measured is therefore smaller, 
while the other class of methods measures the density on a fixed 
scale. 

The environmental parameters used in simulation studies are 
sometimes similar to the ones used for observations, but can also be 
very different. Using a similar definition allows one to directly com- 
pare models and observations. However, with the full (dark mat- 
ter and baryonic) density field available, simulators can also deter- 
mine parameters like the total amount of mass in spheres around the 
galaxy in question. Such quantities might influence the evolution of 
a galaxy, but are difficult or impossible to obtain observationally. 

It is well known that higher mass galaxies preferentially 
live in higher density environments. A coiTelation between halo 
mas s and environmen t al den sity is therefore expected . For exam- 
ple, iKauffmann et alj ( l2004h and ISerrier et all (120 lOl) used semi- 
analytic models of galaxy formation to show how their measure 
of environmental density (number of galaxies within jO,5, I, 2, 3, 
6) /i^'Mpc projected, and a redshift difference less than 1000 km 
s"') correlates with halo mass. They find a good coiTelation with a 
spread of a factor of a few (for small projected cylinders) to a few 
tens (for larger projected cylinders). It is, however, unlikely that 
halo mass is the only characteristic of the en vironment that mat- 
ters. With that in mind, lFakhouri&Mal ( l2009l) tried to construct an 
environmental parameter that does not scale with halo mass. They 
found that the mean over-density in a sphere of 7 Mpc, excluding 
the mass of the halo, gives the most mass-independent parameter of 
the three parameters they studied. They did not quantify the degree 
of correlation, but their plots indicate a weak, but non-negligible 
correlation with host halo mass. Observationally, this quantity can- 
not be determined. As far as we are aware, to date no study has 
found a measure of environment that is independent of halo mass. 



3 POPULAR ENVIRONMENTAL PARAMETERS AND 
THEIR RELATION TO HALO MASS 

In this section we will investigate the relation between some of the 
widely used environmental parameters and the mass of the host 
halo. For the environmental parameters discussed, we will distin- 
guish between the 'ideal case' in which the three-dimensional lo- 
cations and the masses of all galaxies are known (as in simula- 
tions) (Section r3.2l (. and the case in which only projected distances 
and velocity differences can be measured and only luminosities are 
available, as is the case for observations (Section|33). First we will 
briefly summarise the main characteristics of the synthetic galaxy 
populations used. 



3.1 Simulations 

We will compare different environmental parameters using the 
galaxy catalogue c onstructed using the semi-analytic model of 
lOe Lucia & Blaizo3 ( l2007i . see also ICroton et al. 2006h. run o n 
the dark matter-only Millennium Simulation dSpringel et alj2005h . 
The merger trees on which these galaxy models are built, are de- 
rived from Friends-of-Friends haloes (assuming a linking length 
parameter of 0.2) which are subsequently decomposed in bound 
subhaloes by the SubFind algorithm. The Millennium Simula- 
tion follows the evolution of the dark matter distribution using 
2160' particles in a periodic volume of 500 comoving /i ~'Mpc 
down to redshift 0. The model of lOe Lucia & BlaizotI ( l2007h 
uses recipes for the evolution of the baryons inside dark matter 
haloes and is based on the halo merger trees constructed from 
the halo catalogues of the Millennium Simulation. The model 
predicts the galaxies' locations, physical properties such as their 
stellar masses and star formation histories, and observables like 
colours and luminosities. The model is calibrated to reproduce 
the redshift zero luminosit y function in t h e K- a nd i^j -bands. 
iDe Lucia & Blaiz ot (2007 ). iDe Lucia et alj j2007l) . ICroton et al] 
('2006') and lKitzbichler & White! ('2007) showed that this model re- 
produces many other observed properties of the galaxy population 
in the local Universe (e.g. the colour distributions, the stellar mass 
function and the clustering properties). We will only use the z = 
results. We note that the amplitude of fluctuations in the Millen- 
nium Simulation is higher than in the currently favored cosmology 
(trg = 0.9 vs. 0.8). This difference may result in differences in clus- 
tering of dark matter haloes between the Millenniu m Simulation 
and the observed Universe. However, as shown by ICroton et al.l 
(2006), the precursor of the galaxy formation model used here 
reproduced the observed clustering properties very well and the 
model used here does not differ much from its precursor. 

We take into account all galaxies with stellar masses greater 
tha n 10'°/i~'Me. T his is roughly the same lower mass limit 
as iFakhouri & Mai ( [2009) used (they used a total mass of 
1.2xlO'^/i~'M H). This choice is dicta t ed by the resolution limit of 
the simulation. iBovlan-Kolchin et alj i2009i ) showed that the sub- 
halo abundance of haloes in the Millennium Simulation is con- 
verged for subhaloes more massive than about 10"/i"'MQ, roughly 
independent of parent halo mass (as long as the virial mass of the 
parent dark matter halo is larger than 10'^/2''M^). IGuo et alj|201ol) 
also investigated the subhalo abundance convergence of the Mil- 
lennium Simulation and concluded that the halo and subhalo abun- 
dances are conver ged for > 10 '^ 'fe^'Mp. These halo masses 
were matched by iGuo et al. I ( I2OIOI) to the stellar mass function 
from the seventh data release of SDSS from iLi & Whit el (|2009|) . 
from which they conclude that observed galaxies with stellar mass 
M, > lO'^'^/i^'Mo reside in converged haloes. The number of 
neighbours counted in some volume depends on the lower stellar 
mass limit for galaxies in the sample (or, correspondingly, the flux 
limit of the survey), but as we will show, the scalings and correla- 
tions are usually not sensitive to this lower limit. 



3.2 The ideal case: using 3-dimensional distances and masses 

We will use the simplest version of both classes of observationally 
determined parameters: the number of galaxies. A';;, within some 
volume with radius R and the distance to the A'"' nearest neighbour, 
Rn. Parameters derived from these numbers (such as the number 
density of galaxies within that volume) will obey the same conclu- 
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Table 1. Overview of environmental parameters that are frequently used in literature. They are grouped by the dilferent ways of determining out to which 
distance the environment is measured either in observational or simulation studies. The first column specifies the environmental parameter, and the second and 
third column indicate out to what distance the environment is measured and whe the r the minimum mass/lumi nosi ty is fixed or scale s with t he galax y in question. 
The fo urt h column specifies the refe rences for the pap ers: 1: [P ressleJ il98Cll). 2:|Whitmore & Gilmore ( 1 991b. 3: [Gomez et alj i2003h. 4|Post man & Gell3 
il984l ). 5:'Whitmore etal.'('l993'), 6:'Goto et al.' ('2003'), 7 : 'Weinmann et al.' feOOe), 8 : 'Cooper et al. 1 20051) 9: Icooper et al.l j2006l) . 10: [C ooper et al. (200^ 
ll:[B alogh et al. 1 2004a), 12: Balogh et al. l2004b), 13: Baldrv et al. (2006), 14: Bamford et al. (20091) 15:ICassata et al.l j2007h . 1 6:lElli7 on et al. (2010), 17: 
Kovacetal. (2010), 18: Pimbblet et al (2002), 19: Lewis etal. (2002), 20: Blanton et al. (2005), 21: 'Blanton et al.' (■2003b'), 22: Blanton et al. (2003a), 23: 
Hogg et al. (2003), 24: Blanton & Berlind (2007), 25: Berrier et al. (2010), 26: Wilman et al. (2010), 27: Hoga et al. (2004), 28 Kauftmann et al. (2004), 29: 
Park & Cho i (2009), 30: Marker et al. (2006), 31: Hester & Tasitsiomi (2010), 32: Maccio et al. (2007), 33: Maulbetsch et al. (2007), 34: Espino-Brio nes eTal] 
<2007l). 35: Jshivama et al (2008) , 36: Abbas & Sheth ( 2005), 3 7: Grain et al.. 12009.) . 38: .Lemson & Kauftmann, U999.) . 39: .Fakhouri & Ma 120091) . 40: 
iHahn et alj KOOgf) 41: IWang et aljfioOTt) 42: lFaltenbacheil i2010l) 



Parameter 


Distance related parameter value 


Minimum mass/luminosity 


References 


From observations 


(Projected) galaxy number density 


Average of nearest 10 galaxies 


my < 16.5 


1,2,3 






My < -20.4 


3 




Group average 


Mb < -17.5 


4 


Cluster/Group-centric radius 




Mr < -20.5 


5, 6 




- 


Mv < -20.4 


3 






mv < 16.5 


2 




LU llIC Viilai laLllLlA 


r < 17 77 


7 


Projected galaxy number density out 


N = 3, Av = 1000 km s"' 


i? < 24.1 


8, 9, 10 


to the nearest neighbour 


^ = 4,5 


Mr < -20 


11-16 


with a maximum radial velocity 


^ = 4,5, Av= 1000 km S-' 


Mr < -20 


13, 14 


difference Av 


N = 4,5, Av = 1000 km s"' 


Mr < -20.6 


16 




N = 5, Av = 1000 km s"' 


Mr < -20.6 


11 




iV = 5, Av = 1000 km s"' 


Mr < -20 


12 




N = 5, 10, 20, Av = 1000 km s"' 


Iab < 25 


17 




A'= 10 


Mv < -20 


18 




N=m 


/< -24 


15 




N = 10, in clusters 


Mh < -19 


19 


Galaxy number density in sphere 


r ^1 ft^'Mpc 


r< 17.77 


20 


of proper radius r 


r = 8 A^'Mpc, Av < 800 km s"' 


r< 17.77 


21,22, 23 


Number of neighbours in cylinders 


r = 0.1 - 10 /r'Mpc, Av = 1000 km s"' 


Mo.ir - 5Logio/i < -19 


24, 25 


with projected radius r 


r = 0.5, 1, 2 h-'Mpc, Av = 1000 km s"' 


Mr < -20 


26 




r = 1 /i^'Mpc, Av corresponding to 8 Mpc 


r < 17.77 


27 




r = 1 - 10 r'Mpc, Av = 1000 km s"' 


Iab < 25 


17 




r = 2 A^'Mpc, Av = 1000 km s"' 


r < 17.77 


28 


Mass density due to nearest neighbour 


N = I or N for which p is maximal 


Mr.„gb > M,,gal + 0.5 


29 


(p = 3M„gb/4;rr3g^) 


Av =400, 600 km s"' 






Projected galaxy number density in 


10.5,1,2) < R/(h-'Mpc ) < (1,2,3) 


Mr < -20 


26 


annuli 


1 <R/(/i-'Mpc) < 3 


r< 17.77 


28 


From simulations 


Halo mass 




M>2.35xlO"'r'Mo 


30 


Number of neighbours in spheres of radius R 


R = 2r'Mpc 


y„,ax > 120kms-' 


31 


Mass or density in spheres of radius R 


R= 1,2,4,8 /r'Mpc 




32, 33 




i? = 5 r'Mpc 




34, 35 




R = 5, 8 r'Mpc 




36 




R = 7 r'Mpc 




30 




R= 18,25 /r'Mpc 




37 


Matter density in spherical shells 


2 <R/(r'Mpc) < 5 




38, 39, 40 




2<R/(r'Mpc)<7 




30 




RfOF <R< 2 r'Mpc 




30 




^vir < R < 3/?vii- 




41 


Average mass density of suiTounding halos 


^ = 7 


200 < V„,„/kms-'<300 


42 


Distance to nearest halo with minimum mass 




M2/M1 > 3 


35 
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In Fig. [T] we show the correlations between host halo virial 
mass and three definitions of environment: the n umber of galaxies 
withi n 1.5 virial radii (we use the definition of iBrvan & Normanl 
throughout the paper) of the galaxies' host haloes, the num- 
ber of galaxies within 1 ft^'Mpc, and the distance to the fourth near- 
est neighbour (left to right). All three measures of environment are 
strongly correlated with the mass of the host halo. 

If the distance out to which galaxies are counted is scaled to 
the virial radius of the parent halo that the galaxy resides in, then the 
correlation between halo mass and environment is very strong, as is 
shown in the left panel of Fig.[T] Because the region within which 
galaxies are counted grows with halo mass, a more or less constant 
fraction of satellites is counted. A fixed fraction of all satellites is a 
number of satellites that grows roughly linearly with halo mass, re- 
sulting in a very tight correlation . This can be understood in terms 
of the results found bv lGao et al. I( l2004h : the fraction of the mass in 
subhaloes, the distribution of subhaloes and the shape of the sub- 
halo mass function are independent of host halo mass, while the 
normalisation (so the total number of and total mass in subhaloes) 
scales (to first order) linearly with halo mass. The number of sub- 
haloes (and thus satellite galaxies) within a radius that is fixed rel- 
ative to the virial radius therefore grows roughly linearly with halo 
mass. This makes the parameter TV; 5 Rvi^ a very strong measure of 
halo mass. 

A slightly weaker correlation exists between halo virial mass 
and the number of galaxies within a fixed physical distance, as 
shown in the middle panel of Fig.[T](for a distance of 1 /z^'Mpc). 
The upper envelope is populated by the central galaxies in the sam- 
ple, while the satellites form the less tightly correlated cloud above 
the relation for the centrals. 

In the right panel of Fig. [T] we show the correlation between 
host halo mass and the distance to the fourth nearest neighbour, 
(which is very often used observationally, see Table [T). The dis- 
tance decreases with halo mass, because more massive haloes 
are on average found in denser environments. 

The correlation with mass is strongest for ~ l/i"' Mpc, 
which marks the transition from the regime in which the nearest 
4 galaxies are typically part of the same halo {R4 « Ih^' Mpc; 
M > 10'^ /i^'Mq) to the regime in which it resides in another halo 
{R4 » 1/7"' Mpc; M < 10'^/?"'Mo). The transition between the 
two regimes depends on the rank n: for higher ranks, the jump oc- 
curs at higher halo mass. 

The three parameters displayed in Fig.[T|all depend on three- 
dimensional distances. We will now proceed to investigate param- 
eters that can be measured observationally. 

3.3 The realistic case: using projected distances and 
luminosities 

Observationally, we have no access to the three-dimensional sepa- 
rations between galaxies. Instead, we measure distances projected 
on the sky and differences in redshift. Moreover, while luminosi- 
ties are readily available, stellar mass determinations depend on 
SED modelling, which comes with considerable uncertainty. We 
will now investigate to what extent the use of observables weakens 
the correlations compared with the 'ideal cases' discussed in Sec- 
tion [3]2] As is done in many observational studies (see Table [T](, 
we will only make use of galaxies with redshifts that are within 
1000 km s"' of the redshift of the galaxy for which the environ- 
ment is determined. We include both the Hubble flow and peculiar 
velocities in our calculation of the redshifts. For reference, a ve- 
locity difference of 1000 km s"' corresponds to a distance of 10 



/i^'Mpc if the peculiar velocity difference is zero. We will denote 
the parameters using the same symbols as we used for the 3-D dis- 
tance variants, but with lower case letters. For example, r4 denotes 
the projected distance to the fourth nearest neighbour (using only 
galaxies within the redshift difference cut). We only include galax- 
ies with absolute A'-band magnitude smaller than -23, which cor- 
responds to M, « IO'^^/i^'Mq, because in our sample of resolved 
haloes, the galaxy luminosity function shows signs of incomplete- 
ness for fainter galaxies. This results in a slightly smaller sample 
than the one used before. The luminosity function of galaxies with 
M, > lO'^/r'Mo shows signs of incompleteness at magnitudes 
fainter than K = -23. 

Fig. |2] shows the relation between halo mass and the parame- 
ters shown in Fig.[T] but using projected distances and luminosities 
rather than 3-D distances and stellar masses. For all three parame- 
ters the correlations are slightly less strong than in the ideal case. 
Galaxies without any neighbours within the specified distance are 
assigned a number density of log[o(n) = -1. 

Note that the left panel still requires knowledge of the virial 
radius of the host halo of the galaxy and is therefore hard to de- 
termine observationally (we left it in for completeness). The virial 
radius can be estimated if a group catalogue is available, like the 
one bv lYangetailbOOTl) who grouped galaxies using a a friends- 
of-friends like algorithm. The total luminosities of the groups are 
then ranked and matched to a ranked list of halo masses, drawn 
from a halo mass function sampled in a volume equal to that of 
the survey. This procedure results in the assignment of a host halo 
mass to all galaxies in the sample. However, if such a catalogue is 
available, then the halo mass is of course just as well known as the 
virial radius, so using this environmental indicator as a measure of 
halo mass is not very useful. 

In the middle panel of Fig. |2] we show the halo mass as a 
function of the number of galaxies within a projected distance of 1 
/j^'Mpc, with a redshift difference less than ± 1000 km s"' and with 
K < -23. Compared with the 3-D version (middle panel of Fig. [T), 
there are now more low-mass galaxies with a high number of neigh- 
bours. This is due to projection effects. We note that the correlation 
coefficient is still very high (» 0.71), so we can conclude that this 
environmental parameter is a strong indicator of host halo mass. 
The horizontal scatter (scatter in environmental parameter for fixed 
halo mass, so this is not the scatter indicated with the error bars) at 
low halo masses is dominated by projection effects, while at high 
masses the scatter is mainly caused by satellites in the outskirts of 
the halo. The scatter in the environmental indicator is smallest for 
halo masses of about lO''' Ii^'Mq. For a given ni upc/i, the spread in 
halo masses is small for low and high values of the environmental 
indicator (roughly 0.3 dex) and highest for ni Mpc/A ~ 10 (> 0.5 dex 
in halo mass). 

In the right panel of Fig. |2] we show the projected distance to 
the fourth nearest neighbour with K < -23. Because of projection 
effects the bi-modal behaviour visible in the right panel of Fig.[T]has 
been smeared out. The correlation with host halo mass is therefore 
slightly weaker. Because of the discontinuity in the distribution, 
the correlation coefficient is a function of the masses (both galaxy 
stellar mass and host halo mass) of the objects that are taken into 
account. 

3.4 A multi-scale approacii 

IWilman et al] ( I2OI0I) recently measured the number density of 
galaxies in concentric rings of 0.5, 1.0, and 2.0 /j^'Mpc in order 
to investigate trends in the u-r colour distribution of galaxies with 
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Figure 1. The correlations between host halo virial mass and three environmental parameters: the number of galaxies within 1.5 /?vir (left panel), the number 
of galaxies within 1 /i^'Mpc (middle panel) and the distance to the fourth nearest neighbour (right panel) for all galaxies with M, > 10"'/i"'Mo. The numbers 
piinted in the top-left of each panel indicate the Spearman rank correlation coefHcient and the diamond symbols and error bars show the binned medians and 
the 16th and 84th percentiles (for Gaussian distributions this would be ±1(t) of the distribution. All three parameters are measures of halo mass. From left to 
right the correlations with halo mass decrease in strength. 
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Figure 2. As Fig.[T] but now for observable versions of the environmental parameters: We use A'-band luminosities rather than stellar masses for the selection 
of galaxies in the sample and we use projected distances and a cut in redshift difference (of 1000 km s"') rather than 3-D separations for all galaxies with 
K < -23. For the left panel, the virial radii of the host haloes still need to be known. The correlations are slightly weaker than the corresponding ones found 
for the ideal case (Fig. [T), mainly due to projection effects, which make galaxies populate the regions in the plots which were unoccupied in Fig. [T] The 
correlations are, however, still strong. 



environment at several distance scales (for given small-scale den- 
sity, if desired). They included all galaxies from the fifth data re- 
lease of SDSS with magnitude brighter than 17.77 in the r-band 
and with a mean surface brightness within the half-light radius of 
jir < 23.0 mag arcsec"^. The number density of galaxies was de- 
termined in rings with radii fixed in physical coordinates. In this 
approach neither the mass nor the distance out to which the envi- 
ronment is determined scales with the properties of the galaxy in 
question. We therefore expect that these measures of environment 
vary strongly with halo mass. 

We find that the Spearman rank correlation coefficient for the 



density in annuli with halo mass is roughly 0.5, and depends on 
both the width and the radius of the annulus, such that smaller 
radii (< 0.5 Mpc) have larger correlation coefficients and wider 
a nnuli mostly show we aker correlations. The power of the method 
of lWilmanetai] l l2010l) lies in the ability to measure residual trends 
of galaxy properties with large-scale (annular) environment, while 
controlling for the environment on some smaller scale (i.e. the pro- 
jected number density in the inner circle). The samples are con- 
structed by taking all galaxies around which the number density of 
galaxies within the inner radius of the annulus fall within some bin, 
and are therefore comparable to vertical slices through the middle 
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Figure 3. Halo mass as a function of the number of galaxies in annuli with inner and outer radii of 1 and 2 Mpc, respectively. The first panel shows the full 
sample, while the other three panels correspond to three bins in number of galaxies within 1 Mpc, showing the lower, middle and upper one eighth of the total 
distribution from left to right. The galaxies with « = are placed at Logion = -1. Although the range of halo masses in each bin is large, the typical halo 
mass is higher for the bins corresponding to higher values of n \ Mpc/h- The median ;ii Mpc/A for the three bins are 1, 4 and 21, respectively, and the median 
log[Q(Mi,aio/Mo) are 12.3, 12.7 and 14.2. Dividing the sample into narrow bins of 'small-scale environment' generally reduces the strength of the correlation 
between «i_2 Mpc/h and /Whalo> but for large values of ni^^ Mpc/h it actually increases it. 



panel of Fig. |2] From this figure we can see that in such a slice, a 
large range of halo masses is still present. 

As an example, we show in Fig.[3]the correlation between halo 
mass and the number of galaxies in annuli with inner and outer radii 
of 1 and 2 Mpc, respectively, for three narrow bins of the number of 
galaxies within 1 Mpc (projected separation, within a redshift dif- 
ference of 1000 km s"'). Each bin contains 1/8 of all the galaxies, 
where the second panel from the left corresponds to the 1/8 of the 
total galaxy population with the lowest value of n \ Mpc/;i- Similarly, 
the 3rd and 4th panels show the relation between Mhaio and «i_2 Mpc/;i 
for the 1/8 of the galaxies for which ni Mpc/;i is, respectively, in the 
middle and the highest of the values spanned by the total popula- 
tion. The different bins in central number density (i.e. the different 
panels) favour different halo masses, as expected from Fig.|2] Com- 
paring panels 2-4 with the first panel, which shows the results for 
the full population, we see that the correlation between Mhaio and 
wi-2Mpc/;i is much reduced for fixed, low values of «i Mpc/A- How- 
ever, for high values of rii Mpc//i (last panel) the correlation actually 
becomes stronger. 

The trends se en in Fig. [3] ar e a typ ical example of the 'multi- 
scale' approach of IWilman et al.l f201(lt) . Changing the radii of the 
inner and outer edges of the annuli and/or the width of the bins in 
central galaxy number density does not affect the qualitative con- 
clusions drawn from Fig.|3] The correlation of the number of galax- 
ies in annuli with halo mass becomes weaker if very large distances 
from the galaxy in question are taken (5-10 Mpc), but that merely 
reflects the fact that galaxies at such large distances do not have 
much to do with the galaxy in question. 



4 ENVIRONMENT AS A MEASURE OF HALO MASS 

In this section we will study the strength of the correlation between 
several environmental indicators and halo mass in more detail. In 



particular, we will determine which parameter provides the best 
measure of halo mass. 

We expect the correlation between the number of neighbours 
and halo mass to be strongest at some given distance. Taking the 
distance very small will result in strong discreteness effects, while 
taking the distance too large will result in a sample of galaxies that 
does not have much to do with the halo the galaxy resides in. 

In Fig.|4]we show, for two different environmental parameters 
and for two different luminosity cuts, the value of the Spearman 
rank coiTelation coefficient with halo mass, as a function of the 
distance related parameter used to measure the environmental den- 
sity. In the left panel we show the correlation coefficient between 
halo mass and the envirormiental density indicator n, (the number 
of galaxies within a fixed physical distance r projected on the sky 
and within Ai' = ±1000 km s"') as a function of r. An example of 
this type of parameter was shown in the middle panel of Fig.|2]for 
r = i /i^'Mpc. Fig. |4] shows that the correlation first strengthens 
with distance, reaches a maximum at a scale of roughly 1 /i^'Mpc, 
and declines slowly thereafter. The vertical arrows, which indicate 
the median virial radii for the haloes of all galaxies in the sample, 
show that the correlation between iii Mpc/A and Mhaio occurs when r 
is slightly greater than the median virial radius. 

In the right panel of Fig. |4] we plot the Spearman rank 
correlation coefficient between halo mass and environment, now 
parametrised by ^a;, the distance towards the A'* nearest neighbour 
(as in the right panel of Fig. |2] for A' = 4), as a function of the 
rank A'. The correlation coefficients are now mostly negative, as 
a higher density (coiTesponding to a higher halo mass) will result 
in a smaller distance towards the A'"' nearest neighbour. However, 
if we only consider galaxies with K < -25.5, then halo mass is 
an increasing function of the distance to the nearest neighbour as 
the neighbour needs to be outside the galaxy itself, and more mas- 
sive galaxies tend to be larger. Taking more neighbours, i.e. higher 
values of A', gives an anti-correlation that, for high-mass galax- 
ies, becomes stronger for larger numbers of neighbours. For lower 
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Figure 4. The strength of the correlations between halo mass and two of the environmental indicators that can be used straightforwardly for observations for 
two samples of galaxies with the luminosity Hmits indicated in the legends. In the left panel we plot the Speannan rank correlation coefficient between halo 
mass and the number of galaxies within a given projected physical distance r (and with a cut in redshift difference of 1000 km s"') as a function of r. The 
arrows show the value of the median virial radius of the haloes of all galaxies in the sample with the corresponding colours. The right panel shows the Spearman 
rank coiTelation coefficient between halo mass and the projected distance to the nearest neighbour as a function of the rank N. The correlation coefficient 
is typically negative, because more massive galaxies have their A?"" nearest neighbour closer by. The arrows indicate the median number of neighbours within 
the virial radius of the haloes above the indicated flux limit. The parameters n,- with r > 1 Mpc and with N ~ iO are very good measures of halo mass for a 
wide range of luminosities. 



mass galaxies (K < -23) the correlation between and Mhaio is 
strongest for N x 3-4, but does not weaken much for larger values. 
The median number of neighbours within the virial radius, above 
the same luminosity cut is indicated with the arrows. 

We conclude that n,. and r^. are both good measures of host 
halo mass, provided that n, is measured at r > r^i^ and that A' is suf- 
ficiently large (A' x 10). If the host halo mass, and thus the virial 
radius, is not known a priori, then it is better to take r larger (r > 1 
/j^'Mpc), as the correlation rapidly weakens towards smaller dis- 
tances and declines only slowly with increasing distance. Although 
the exact values of the correlation coefficients depend on the lower 
limits in halo mass and galaxy luminosity used in the study, the 
qualitative result that the scale for which the correlation is strongest 
is comparable to the median virial radius, is robust. The qualitative 
results are also expected to be similar at higher redshift. 

Observational studies often contrast stellar mass and environ- 
ment. A dependence on environment for fixed stellar mass is as- 
sumed to reflect nurture rather than nature. However, Fig. |5]shows 
that breaking up the sample in small bins of /f -band magnitude for 
the galaxy under consideration (in the neighbour search we still in- 
clude all galaxies with K < -23), which correlates well with stellar 
mass for central galaxies, does not reduce the strength of the corre- 
lations between environment and halo mass. Thus, effects attributed 
to environment, even in studies that fixed stellar mass, are mostly 
due to halo mass. This is important, as halo mass, at least for cen- 
trals, is considered to reflect nature rather than nurture in models of 
galaxy formation. 

As we will show below, using A^-band luminosity as a proxy 
for (virial) mass works well. Guided by the left panel of Fig. |2] one 
might expect that we can improve on as a measure of halo mass 
if r scales with lJ/' . We have tried this, but the coiTelation between 
halo mass and environment does not get stronger (or it gets slightly 
weaker, with correlation coefficients of 0.65 - 0.7). In the range of 



halo masses for which we could test it (any range between 10'^ and 
IO'^'^A^'Mo) the correlation is stronger if a projected distance of 
1 Mpc is used than if r oc lJ/' is used. Specifically, we tried r = 
l/i-'Mpc (Ljf/Lo)"-\ with Lq = loHO-^.n 0.11.5.12.01^^ We therefore 
conclude that using a fixed physical projected distance is safe, and 
easier in practice than a distance scaling with luminosity. We thus 
advise to use n,. with r of the order of r > Ry,y, if a measure of 
halo mass is desired. For most observed samples of galaxies r ~ I 
Mpc will do, but by iteration better values can be obtained: use 
r = 1 /z^'Mpc, calculate the halo virial radii from the environmental 
indicator (using the parametrisation given in Appendix |A) and then 
iterate if the virial radii strongly deviate from 1 Mpc. 

In Appendix |A] we provide polynomial fits for the halo mass 
as a function of several environmental parameters for several flux 
limits, which can be used to obtain halo masses from observed sam- 
ples of galaxies with measured environmental indicators. Note that 
the fits are for z = and may differ for higher redshift samples. 



5 ENVIRONMENT INDEPENDENT OF HALO MASS 

In the previous section we have considered which environmental 
parameters are the best measures of halo mass. In this section we 
will construct new measures of environment, both for use with sim- 
ulations ( 95. It and with observations ( 95.21 . that are highly insen- 
sitive of halo mass. 



5.1 Halo mass independent parameters for simulations 

All the environmental parameters we have looked at so far corre- 
late with halo mass. The lower mass/luminosity limit of galaxies in- 
cluded as possible neighbours was set equal to the resolution limit 
of the simulations, or the flux limit of a survey. As we saw in the 
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Figure 5. The same as Fig.|4] but now for three bins in absolute S'-band magnitude for the galaxy under consideration (for the neighbour search all galaxies 
with K < -23 are taken into account). The numbers in between the brackets indicate the number of galaxies in the sample. The correlations between the 
environmental parameters and halo mass are insensitive to the A'-band luminosity (and thus to stellar mass) except for small values of r and N. Even for narrow 
bins in K, the environmental parameters correlate strongly with halo mass. 



left panels of Figs. [T] [J] andU the correlation is strongest, and al- 
most linear with halo mass, if we count neighbouring galaxies out 
to the virial radius of the host halo of the galaxy in question. Per 
unit halo mass, this galaxy number density (either projected or in a 
spherical region) is therefore roughly constant. This also holds for 
dark matter subh aloes in high-resolution simulations, as shown by 
lGaoetal.l l l2004h . 

In order to obtain an environmental indicator that is indepen- 
dent of halo mass, we expect that we will have to scale out both the 
mass/luminosity of the galaxy and the length scale in question. We 
define D^j to be the three-dimensional distance to the A'* nearest 
neighbour with a viral mass that is at least / times that of the halo 
under consideration, divided by the virial radius of the A'"' nearest 
neighbour: 



D, 



N.f 



(1) 



vir. ngb 



where the subscripts 'ngb' and 'halo' indicate, respectively, the 
neighbour of the halo under consideration and the halo itself. Ob- 
serve that the adjustable parameters A' and /, as well as Djv,/ itself, 
are all dimensionless. Because the tidal force due to the A'* nearest 
neighbour scales as M„^f,/Rjj oc (./?vii.ngb/^«)^> the parameter D^, ; 
scales with the tidal force to the power -1/3. This makes D^j a 
natural envirormiental parameter with a clear physical interpreta- 
tion.^ 

IPark & Choil (12009') found in an observational study that the 
physical properties of galaxies are sensitive to the distance to the 
nearest neighbour, normalised by the virial radius of that neighbour 
They estimated virial properties of the galaxies from the luminosi- 
ties. We will use a similar approach in the next section. 

The colour scale of Fig. |6] shows the distribution of haloes at 
z = in the Di i - Mhaio plane. The curve shows the median Di i 
in bins of halo mass. Although the median Di f in the sample is 
different for different /, halo mass is always independent of D, ir- 
respective of the factor /. The weak correlation that appears at the 
highest values of Di j and Mt^.^i^ is caused by the fact that these 
massive haloes are on the exponential tail of the Schechter-like 
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Figure 6. Halo mass as a function of the environmental parameter D[ i 
(eq.[T). The colour scale gives the distribution for all central galaxies in the 
sample, while the soHd line is the median halo mass in bins of Z)| i . The 
median relation is very flat. The correlation coefficient of this parameter 
with halo mass is 0.07 (for congelation coefficients as a function of rank, 
see Fig.|2). Except for the very high Di i end, where the median halo mass 
is sufficiently high for the haloes to be on the exponential tail of the mass 
function, Di [ and halo mass are uncorrelated. 



halo mass function. Large-scale structure is no longer self-similar 
in that regime, causing a slight positive correlation between D^j 
and halo mass. We have verified (by inverting the axes) that for 
masses M ^ M, (where M, is the mass at which the Schechter- 
like halo mass function transits from a power law into an exponen- 
tial fall-off), where the mass function is a power law (and therefore 
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Figure 7. The Spearman rank correlation coefficient between halo mass and 
the environmental indicator Da/,/- (see Eq.[T) as a function of the rank N, for 
/ = (1/10, 1, 10). Higher values for / and N result in stronger correlations. 
If haloes can be reliably identified for masses lower than the lowest mass 
one wants to know the environment for, then using a low value of / (e.g. 
/ = 0.1) gives negligible correlations with mass. For / = 1 the environ- 
mental parameter will be measurable for all haloes that can be identified 
and this value of / still gives only very weak correlations with halo mass, 
particularly for low values of N. 



scale free), the correlation is very weak. For higher masses, the ex- 
ponential cut-ofF of the Schechter-like halo mass function imposes 
a mass scale. For values above roughly /"'M,, the insensitivity to 
mass breaks down and a weak positive correlation between halo 
mass and Df^j appears. 

In Fig. |7] we show the correlation coefficients between halo 
mass and Df^ f as a function of the rank N for three different val- 
ues of the mass ratio of galaxies counted as neighbours and the 
mass of the halo under consideration, / = jl/10, l,10j. For all val- 
ues of / the correlation between the rank A' and host halo mass 
increases with the ra If an environmental indicator is desired 
that is insensitive to halo mass, N = 1 is therefore a good choice. 
The correlation is weaker for lower values of /. For / < 1 the 
environmental indicator cannot be determined for the full resolved 
sample of haloes (as halo masses need to be at least M > f^^M^^^, 
with Mpes the resolution limit, in order to resolve all possible neigh- 
bours). We therefore advise to take / = 1, as then the parameter 
can be defined for all galaxies in the sample and it still gives only 
a very weak correlation with halo mass. If in a sample of haloes 
some of the studied properties demand a much more stringent res- 
olution limit (e.g. if detailed halo profiles need to be fitted), and 
if haloes of much lower mass are resolved in terms of their virial 
mass and position, then one should use values of / < 1 , e.g. 0. 1 , for 
which the correlation between halo mass and environment becomes 
vanishingly small. 

If, in the definition of D/yj, we replace the virial radius of 
the neighbour by the virial radius of the halo under consideration 



' Note that the correlation does not necessarily vanish for A' — > oo, because 
neighbours must be at least / times as massive as the halo under considera- 
tion. However, for very large N the parameter D^.f no longer characterises 
the environment of the galaxy as the distance to the nearest neighbour 
diverges with N. 



(thereby losing the connection to the tidal force due to the neigh- 
bour), the correlation between halo mass and environment becomes 
even slightly weaker (e.g. a Spearman rank correlation coefficient 
of 0.04 instead of 0.07 between halo mass and fi.i). As using the 
virial radius of the neighbour gives a more intuitive external envi- 
ronmental parameter, we prefer to use the virial radius of the neigh- 
bour. 

We conclude that the parameter D/^ f, with N = I and / < 1, 
results in an intuitive environmental parameter that is very insensi- 
tive to halo mass. We do note, however, that in order to calculate 
this halo mass independent environmental indicator, one needs a 
measure of the virial mass of the host halo. From simulations these 
can easily be obtained. For observed samples of galaxies virial 
masses can be estimated using the environmental indicators that do 
correlate strongly with halo mass, as described in the previous sec- 
tion and detailed in AppendixlAl In the next section we will present 
an environmental indicator that can be obtained directly from ob- 
servations and that is also insensitive to halo mass. 



5.2 Halo mass independent parameters for observed samples 
of galaxies 

The environmental parameter Dj^ f defined in the previous section 
depends on the masses of the haloes, which are generally not read- 
ily available for observed samples of galaxi es. Halo masses can 
e.g. be estimated by abundance matching (e. g.'Kravtsov et al.'2004|; 
I Vale & OstrikeJ |2004 IConrov et al] I2OO 6: Conroy & Wechsler| 
'2009^, where a sorted list of galaxy luminosities is matched to a 
sorted list of halo masses (from either an analytic halo mass func- 
tion or a simulation), with or without scatter in the one-to-one re- 
lation. However, this technique does not work for satellite galaxies 
and it is generally not known which galaxies are centrals and which 
are satellites. If a group catalogue is available, then the results can 
be improved by summing the luminosities of all galaxies in a group 
and matching t he sorted list of g roup luminosities to the sorte d list 
of ha l o masses Yang et al. (e.g. l2003h : Ivan den Bosch et (e.g. 
|2003|); lYang et al.l (e.g. \20m . Finally, as we have shown in Sec- 
tion |4] it is possible to use common environmental parameters to 
estimate halo mass. We provide detailed instructions for doing so 
in the Appendix. 

Since it is easier to work with environmental parameters that 
do not require knowledge of halo masses and because Df^j- also re- 
quired knowledge of 3-dimensional distances, we set out to formu- 
late an environmental parameter that can be easily determined ob- 
servationally and that is as insensitive of halo mass as possible. We 
let the definition of D^j- guide us. We know that we have to scale 
the minimum masses/luminosities of the galaxies that are taken into 
consideration in the search for neighbours to be a fixed fraction of 
the mass/luminosity of the galaxy under consideration and that we 
have to scale the distance to the neighbours to some typical length 
scale associated with the neighbour. 

We use an observable, the A'-band luminosity Lk, instead of 
mass. Luminosity is easier to measure and does not require the 
modeling of the spectral energy distribution of the galaxy. We use 
the A'-band because in the very red optical bands and in the near- 
IR the correlation between luminosity and stellar mass is strongest 
(aside from the uncertainties arising from the treatme nt of ther- 
mally pulsing asymptoti c giant branch stars, see e.g. iMarastonI 
l200alToninietal.ll201(]|) . Instead of 3-dimensional distances we 
will use projected distances, considering only neighbours with red- 
shifts that are within 1000 km s"' of the galaxy for which we are 
measuring the environment. Instead of dividing by the virial ra- 
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Figure 8. The host halo masses of galaxies versus the dimensionless environmental parameter di o (see Eq.|2j, which is easy to measure from observations. 
In the left panel we show all galaxies in the sample, for which the Spearman rank coiTelation coefficient is -0.28. The middle and right panels show the same 
relation for central and satellite galaxies, respectively. For these sub-samples the Spearman rank correlation coefficients are 0.09 and -0.35, respectively. The 
correlation between M[,jiio and di o is weak and nearly vanishes for centrals. 



dius of the neighbour, we divide by L^, as i?vii M''^ and we 
expect that will scale roughly as Mi,aio, at least for central galax- 
ies. Note that this last assumption is known to break down (e.g. 
iKravtsov et al. 2004; Coorav 2006; Conrov & Wechsleill2009l) , but 
we find that using more complicated functional forms to fit the re- 
lation between i^vir and {Lie) does not significantly change the cor- 
relation between environment and halo mass. 

Our environmental indicator df^,,, then becomes 

"■'""0.58/i-iMpcAl.4xlO"Lo/ ^' 

where the subscript 'ngb' again denotes the neighbour of the galaxy 
in question, m is the difference in magnitudes (corresponding to a 
ratio in luminosity/mass, a positive m implies that the neighbours 
must be brighter) between the galaxy in question and the galax- 
ies counted as possible neighbours, and K is the absolute K-band 
magnitude. Although the coefficients do not affect the correlation 
with halo mass, we have divided the projected distance by /?vir.i3 = 
0.58r' Mpc(Lk/1.4 X 10" Lq)""' rather than just by lJ/', where 
^vir,i3 = 0.58/2"'Mpc is the projected virial radius of the 'reference 
mass' of lO'^/j^'lVIo, to ensure that the actual values of rfiv,m re- 
tain some intuitive, physical meaning. If /?vir,i3(i'A'/l-4 x 10"Lo)''' 
were the virial radius, then the external environmental indicator 
df^„, could be described as the projected distance to the A'"' nearest 
neighbour that is at least m magnitudes brighter than the galaxy we 
are measuring the environment of, normalised to the neighbour's 
projected virial radius. 

The colour scale in the left panel of Fig.[8]shows the distribu- 
tion of galaxies in the Mj^io - di o plane. We include all galaxies in 
the catalogue with K < -23. The Spearman rank correlation coeffi- 
cient between halo mass and environment is -0.28, which indicates 
a weak anti-correlation. 

The middle and right panels of Fig. [8] show the distribu- 
tion of central galaxies and satellites, respectively. For these sub- 
samples the Spearman rank correlation coefficient between di o and 
halo mass are 0.09 and -0.35, respectively. Thus, the weak anti- 
correlation seen in the left panel is mostly due to the inclusion of 
satellites. IVIore massive haloes contain more satellites above a fixed 
magnitude limit. Hence, for satellites in more massive haloes, the 
distance to the nearest brighter galaxy, which can be another satel- 
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Figure 9. The Spearman rank congelation coefficient between halo mass 
and cInjii as a function of the rank N, for m = {-2.5,0,2.5) magnitudes 
(i.e. neighbours with A'-band luminosities that are at least 0.1, 1.0, and 10 
times as large as that of the galaxy under consideration). If an environmental 
indicator that is insensitive to halo mass is desired, then rfi 2.5 and d^o for 
any N < 10 are all good choices. 



lite or the central galaxy of the same halo, is typically smaller rela- 
tive to the virial radius of the halo. 

The parameter shown in Fig.[8]requires the difference in red- 
shift between the neighbour and the galaxy under consideration to 
be smaller than 1000 km s"' . Without this cut in redshift difference 
the correlations become stronger. Taking into account only galax- 
ies within a redshift window is important, but the result is not very 
sensitive to the precise velocity cut. 

The dependence of the correlation between host halo mass and 
dN,m on the rank A' is shown in Fig.|9l for three different values of 
m. We have chosen to show m = j-2.5, 0, 2.5) magnitudes, because 
a magnitude difference of 2.5 corresponds to a luminosity ratio of 
10, similar to the mass ratio of 10 that we used earlier. If neighbours 
that are less luminous than the galaxy under consideration are al- 
lowed to count as neighbours (i.e. m negative), then the sample for 
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which the djvm can be determined is smaller than the total sam- 
ple of galaxies (because all possible, lower mass neighbours need 
to be resolved) and the typical haloes the galaxies are in are more 
massive. For m = -2.5, i.e. neighbours that are at least 0.1 times 
as bright, we find a weak anti-correlation that slowly decreases in 
strength with increasing N, from 5 ~ -0.2 for A' = 1 to 5 ~ -0.1 
for A' = 10. For m = the anti-correlation is weakest for A' = 1, 
the case shown in Fig. [8] For in = 2.5, i.e. neighbours that are at 
least 10 times more luminous than the galaxy under consideration, 
we find a very small, positive correlation for A' = 1 (5 ~ 0.07), 
which turns into a anti-coiTelation for N = 2 that quickly increases 
in strength with increasing A'. Hence, if an environmental indicator 
that is insensitive to halo mass is desired, then di j.s and d/^o for any 
A' < 10 are all good choices. 



6 CONCLUSIONS 

The properties of observed galaxies and dark matter haloes in sim- 
ulations depend on their environment. The term "environment" has, 
however, been used to describe a wide variety of measures that may 
or may not correlate with each other. Useful measures of environ- 
ment include, for example, the distance to the A'"' nearest neigh- 
bour, the number density of objects within some distance, or, for 
the case of galaxies, the mass of the host dark matter halo. In this 
paper we carried out a detailed investigation of several environmen- 
tal parameters which are popular in the (observational) literature, 
focusing in particular on their relationship with halo mass. 

We measured the environmental indicators from the syn- 
th etic galaxy c atalogues produced using the semi-analytic model 
by De Lucia & B laizot (2007), built on the Millennium Simulation 
jSpringel et al.ll2005l) . This model reproduces the number density 
and clustering properties of observed galaxies in the low-redshift 
Universe. All our results hold for galaxies with a stellar mass in ex- 



cess of 10'"/2"'Mg in haloes above 10'- A^'Mq at z ~ 0. At lower 
masses and higher redshifts the results are expected to be qualita- 
tively similar. 

We showed that it is of crucial importance to realise that the 
degree to which environmental parameters measure host dark mat- 
ter halo mass, is determined by (1) whether the distance out to 
which the environment is measured scales with some typical length 
scale (e.g. the virial radius of the halo hosting the neighbour) and 
(2) whether or not the minimum mass/luminosity that the neigh- 
bours are required to have is fixed in absolute terms or relative to the 
mass/luminosity of the galaxy in question. Specifically, we found 
that 

(i) All frequently used environmental indicators (i.e. some func- 
tion of the distance to the A'"' nearest neighbour or the number of 
galaxies within some given distance, either using three-dimensional 
distances or using projected distances for all galaxies within some 
radial velocity difference) correlate strongly with halo mass. This 
remains true if only galaxies within some narrow range of (A"-band) 
luminosities are considered. 

(ii) For the number of galaxies within a distance r, iir, the cor- 
relation with halo mass is strongest if we set r equal to 1.5-2 virial 
radii. The virial radius is generally difficult to measure from obser- 
vations and knowing it would remove the need to measure halo 
mass, but the correlation with halo mass is nearly as strong for 
galaxy counts within ~ 1 Mpc. 

(iii) The strength of the anti-correlation between the distance to 
the A'* nearest neighbour, ta-, and halo mass is nearly constant for 
N ^ 2 and only slightly weaker for A' = 1. The relation between 



halo mass and is slightly weaker than that between halo mass 
and n,. if r is taken to be similar to the virial radius. 

(iv) Both and correlate more strongly with halo mass if the 
neighbours are required to be more luminous or massive. 

We have shown that it is possible to construct environmental 
parameters that are highly insensitive to halo mass by using only 
dimensionless quantities. For the case of dark matter haloes in nu- 
merical simulations this can for example be achieved by measuring 
Dnj, the distance to the A'* nearest halo, that is at least / times 
as massive as the halo under consideration, divided by the virial 
radius of that neighbour. Dividing by the virial radius of the halo 
hosting the galaxy itself yields even slightly weaker correlations 
with halo mass, but dividing by the virial radius of the neighbour 
gives the indicator an intuitive interpretation: the tidal force due 
to the neighbour scales as l/D^^ . The correlation with halo mass 
becomes weaker if the minimum mass required for neighbours is 
lower (i.e. for lower /). These environmental parameters are, how- 
ever, only insensitive to halo mass for haloes that are not on the 
exponential tail of the mass function. 

In the case of observations, we usually only know a position 
on the sky, some rough indication of the distance along the line of 
sight, and the flux or luminosity in some waveband. We showed 
that analogous environmental measures that are highly insensitive 
to halo mass can also be constructed using only the A'-band lu- 
minosities, projected distances on the sky, and a maximum radial 
velocity difference for neighbours. Specifically, the parameter di^„,, 
defined as the projected distance to the A'"' nearest galaxy that is at 
least m magnitudes brighter and that is within a radial velocity dif- 
ference of 1000 km s"', divided by the AT-band luminosity of that 
neighbour to the power one third, correlates only very weakly with 
host halo mass for suitable choices of A' and m (e.g. m = 2.5 and 
A' = 1 or m = and any A' < 10). The correlation vanishes nearly 
completely for samples that only contain central galaxies. 

In summary, when measuring environments for (virtual) ob- 
servations, we advise to make use of both a halo mass independent 
measure and a measure that is highly sensitive to halo mass. For 
purely theoretical studies the halo mass is already known and we 
therefore advise to use an environmental parameter that is insensi- 
tive to halo mass. The following parameters are good choices: 

• Insensitive to halo mass; for simulations: The distance to the 
nearest (main) halo that is at least / times more massive than the 
halo in question, divided by the virial radius of that neighbour. The 
choice / = 1 works well, but if resolution permits it, smaller values 
yield even weaker correlations with halo mass. Dividing instead by 
the virial radius of the halo itself gives a slightly weaker correlation 
with halo mass, at the expense of losing the intuitive definition in 
which the environment relates to the tidal force due to the neigh- 
bour. 

• Insensitive to halo mass; for obsen'ations: The parameter di o, 
as given by Eq.|2] The correlation with halo mass is weaker if satel- 
lites are excluded. 

• Sensitive to halo mass; for observations; The number of 
brighter galaxies within a projected distance of ~ 1 /i^'Mpc, within 
a redshift window corresponding to Av < 1000 km s"'(ni Mpc/h)- 
Even better would be to subsequently iterate the following two 
steps until the procedure converges: (i) check what the correspond- 
ing halo masses are using the relations between n,. and halo mass 
given in AppendixlA] (ii) adapt the maximum projected distance to 
1.5 times the typical virial radius of the haloes in the sample. 

Many studies have measured galaxy properties as a function 
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of both stellar mass and environment. We have shown that the en- 
vironmental indicators used by most authors are effectively mea- 
sures of halo mass, even for fixed A^-band luminosity, which is a 
proxy for stellar mass. While halo mass is a perfectly valid measure 
of environment, and may be particularly relevant for satellites, we 
note that because stellar mass is also expected to correlate strongly 
with halo mass, these studies may not have separated "internal" 
and "external" influences as well as one might naively think. The 
work presented here will enable future observational and theoreti- 
cal studies to disentangle the effects of halo mass (internal environ- 
ment) from those of the external environment. This may eventually 
tell us whether halo mass is the only important driver of the physics 
governing galaxy evolution. 
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APPENDIX A: OBTAINING THE HALO MASS FROM 
ENVIRONMENTAL PARAMETERS 

In this Appendix we describe the tables we provide in 
the electronic edition of the paper and on the website 
http://environment.marcelhaas.com These tables provide fitting 
functions for the halo mass as a function of different environmental 
indicators. 

We will use environmental parameters that can be obtained 
directly from observations and those that need an iterative scheme 
outlined below, in order to estimate the virial radius of the host 
haloes of the galaxies in question. We provide the parameters cor- 
responding to third order polynomial fits for the halo mass as a 
function of the environmental indicators. We fit a function of the 
form 

logio[Mh,,o(r 'Mo)] =A + BP + CP- + DP^ (Al) 

Where P indicates the logarithm (base 10) of the environmental 
parameter in question. For all indicators we fit to the median halo 
masses in bins separated by IS.P = 0.25. 

The fitted values for the four polynomial coefficients 



(A, B, C, D) are given in the online tables (see sample table lAlt 
for nine different environmental parameters (ri, r^, rio, no.5Mpc/A> 

"I Mpc/A, «2 Mpc/A, , «l Rvir, "1.5 Rvir and «2 Rvir)- 

Because the distribution in halo masses at fixed environment 
is neither Gaussian nor symmetric, we provide similar fits to the 
16* and 84"" percentiles of the data set (these percentiles would 
correspond to ±\cr for Gaussian distributions), such that at all val- 
ues for the environment the median and 16* and 84* percentiles of 
the distribution of halo masses can be found. 

We provide tables for K-band magnitude limited samples with 
galaxies brighter than a: = {-23,-23.5,-24,-24.5,-25,-25.5) 
and tables for stellar mass limited samples with 
log,JM.(/;-'M,™)] > 110,10.2,10.4,10.6,10.8,11}. In the 
neighbour search we include all galaxies above the given flux or 
stellar mass limit, with a maximum redshift difference of 1000 
km/s between the neighbour and the galaxy we determine the 
environmental parameter for. 

For all parameters and (stellar mass or flux limited) samples 
we provide fits for the halo mass including all galaxies, or in small 
bins of stellar mass or A'-band magnitude. Note that the obtained 
halo masses may become unreliable below lO'^Msi,,,, the conver- 
gence limit of the simulation, and that the fits were derived for 
z = Q. For higher redshifts and lower halo masses, a similar analysis 
would need to be carried out. 

Al Iterating towards the best halo mass estimator 

As we have shown in Section [T2l the strongest correlation between 
halo mass and environment is obtained when galaxies are counted 
within a distance that scales with the virial radius of the halo host- 
ing the galaxy under consideration. In order to do so, an estimate 
of the host halo mass is required, which can be obtained from the 
relations described earlier in this Appendix. Using 

= 0.27 /r'Mpc f ] ' — , (A2) 

which is the relation that was used in this work to obtain virial radii, 
an estimate for the virial radius can then be obtained. Here z is the 
redshift, which is zero throughout this work. 

The best estimate of the halo mass can then be found by mea- 
suring the projected number of neighbours within a given multi- 
ple of the virial radius (with the same cut in radial velocity differ- 
ence), as shown in Section 1331 In the tables we provide the same 
third order polynomial fits, but for the relation between halo mass 
and «i Rvir> "1.5 Rvir ^u )i2Rviii Well as the corresponding (higher) 
Spearman rank correlation coefficients. 

This procedure for obtaining a better estimate for the halo 
mass can then be iterated towards a reliable estimate for the halo 
mass, including the spread in halo masses at fixed environment 
(note that this spread is very small for high-mass haloes if the 
neighbours are counted within a multiple of the virial radius of 
about one.) 

Finally, we caution that these halo masses are measured from 
the Millennium Simulation, which uses the WMAP first-year re- 
sults for the cosmology, which has (among other differences) too 
large an amplitude of fluctuations (erg = 0.9). This means that for 
a given galaxy luminosity, the haloes will be slightly too massive. 
How this affects the relations between environment and halo mass 
is not clear. 
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Table Al. A sample table from the complete tables provided in the online edition and on our website fhttp://environment.marcelhaas.com't. This table contains 
the fits for halo mass as a function of environment for all galaxies brighter than K = —23 and is shown here for two of the environmental parameters. The first 
two columns give the upper and lower flux limits for which the fits are made {K = -999 indicates no upper flux limit was used). Note that in the neighbour 
search all galaxies brighter than K = -23 are used (within a redshift difference of 1000 km s"'). Columns 3-6 are the coefficients A - D of Eq. lAll for the 
median relation (subscript 'med'), and columns 7-10 and 11 - 14 are the coefficients A - D of EQ. lATI for the 16"" and 84* percentiles (subscripts '16' and 
'84'), respectively. The last column gives the Speannan rank correlation coefficient between host halo mass and the environmental parameter Note that all 
lines for a given environmental parameter, except for the first and last, can be used to determine the host halo mass from the environment for fixed K band 
luminosity (and for fixed stellar mass in the corresponding tables). 
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0.55 


0.37 


12.96 


-1.5 


0.5 


0.6 


-0.24 


-23.5 




-24.0 


12.4 


-0.55 


0.57 


0.38 


12.01 


-0.20 


0.55 


0.29 


12.97 


-1.1 


0.5 


0.5 


-0.48 


-24.0 




-24.5 


12.2 


-0.65 


0.90 


0.57 


11.89 


-0.24 


0.87 


0.45 


12.73 


-1.3 


0.9 


0.7 


-0.54 


-24.5 




-25.0 


12.2 


-0.66 


0.79 


0.48 


11.89 


-0.31 


0.77 


0.41 


12.70 


-1.1 


0.7 


0.6 


-0.58 


-25.0 




-25.5 


12.3 


-0.67 


0.51 


0.30 


12.09 


-0.33 


0.34 


0.16 


12.73 


-0.9 


0.5 


0.4 


-0.65 


-25.5 




-999.0 


12.4 


-0.88 


0.75 


0.42 


12.17 


-0.46 


0.61 


0.24 


12.69 


-1.2 


0.7 


0.5 


-0.70 


P = logio[r4(/i" 


'Mpc)] 






























-23.0 




-999.0 


12.8 


-0.87 


0.51 


0.30 


12.37 


-0.89 


0.49 


0.33 


13.40 


-0.9 


0.4 


0.3 


-0.70 


-23.0 




-23.5 


12.8 


-1.02 


0.56 


0.42 


12.34 


-1.01 


0.55 


0.41 


13.38 


-1.1 


0.5 


0.4 


-0.72 


-23.5 




-24.0 


12.8 


-0.85 


0.47 


0.29 


12.38 


-0.92 


0.47 


0.37 


13.44 


-0.9 


0.4 


0.3 


-0.70 


-24.0 




-24.5 


12.8 


-0.79 


0.57 


0.29 


12.29 


-0.88 


0.58 


0.36 


13.37 


-0.9 


0.5 


0.3 


-0.65 


-24.5 




-25.0 


12.6 


-1.16 


0.73 


0.49 


12.29 


-0.97 


0.60 


0.43 


13.21 


-1.1 


0.6 


0.5 


-0.62 


-25.0 




-25.5 


12.7 


-1.08 


0.50 


0.35 


12.33 


-0.90 


0.58 


0.31 


13.18 


-1.3 


0.4 


0.4 


-0.68 


-25.5 




-999.0 


12.9 


-1.43 


0.31 


0.47 


12.58 


-1.26 


0.47 


0.48 


13.36 


-1.4 


0.1 


0.3 


-0.83 



